Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



Realization of Strong Coupling Fixed Point in Multilevel Kondo Models 
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Impurity four- and six-level Kondo model, in which an ion is tunneling among four- and six- 
stable points and interacting with surrounding conduction electrons, are investigated by using 
the perturbative and numerical renormalization group methods. It is shown that purely orbital 
Kondo effects occur at low temperatures in these systems which are direct generalizations of 
the Kondo effect in the so-called two-level system. This result offers a good explanation for the 
enhanced and magnetically robust Sommerfeld coefficient observed in SmOs4Sbi2 and some 
other filled-skutterudites. 
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Since the two-level Kondo effect was studied about 
three decades ago,^'^ the physics arising from the off- 
center degrees of freedom of ions have attracted much at- 
tention. Recently, off-center motions of ions in clathrate 
compounds^ and filled skutterudites^ have been iden- 
tified by ultrasonics. One of the recent remarkable ex- 
perimental results is a realization of magnetically robust 
heavy fermion in SmOs4Sbi2, in which the linear temper- 
ature coefficient of the specific heat at low temperatures 
amounts to 820mJ/(mol-K2) even up to 8-15 T.^'"^ Since 
the Sm ion is considered to be in a mixed valence state, it 
is mysterious why such a heavy effective mass is realized. 

Characteristic properties of these materials are the ex- 
istence of thermally excited rattling modes if an ion in 
the "cage". In the case of PrOs4Sbi2, it is claimed that 
there exist degenerate Fj^ (cubic Oh) modes of the off- 
center motion of the ion (Pr'^+) to explain the ultra- 
sonic dispersion in Cn — C12. (This point has been still 
controversial. There are two different results in (744^'^). 
However, textbooks for quantum mechanics tell us that a 
ground state is a singlet with a symmetric wavefunction 
and there exists no node in a spinless one-body prob- 
lem. From this view point, the interpretation above is 
not easily accepted. 

In this paper, we investigate an impurity four- and six- 
level Kondo model, in which an ion is tunneling among 
four- and six-stable points. We take into account inter- 
actions between the ion and the surrounding conduc- 
tion electrons, and use the perturbative renormalization 
group (PRG) method to see low energy properties of the 
system. Since we want to investigate the effects arising 
from the non-magnetic origin, we assume the conduction 
electron to be spinless. We restrict ourselves mainly to 
the six-level model (see ref. 8 for the detailed PRG and 
numerical renormalization group (NRG) discussions of 
four- level model). We will make a brief comment about 
the results of the four- level model in the last part. 

In order to make the model simple, we consider a Gaus- 
sian wavefunction of the ion (jJf (x) at each off-center po- 
sition i. 
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where a and are the width of the Gaus- 

sian and coordinate of off-center positions x^ — 
(±a,0,0), (0, ±a,0) and (0,0, ±a), respectively (a and 
a are measured in the unit of inverse Fermi wavenumber 
fc^^ ). Using 0°(x), we can evaluate the transfer inte- 
gral S directly (we keep only nearest neighbor terms). 
Since S = exp(— a^/(2a^)) ^ 0, the energy eigenstates 
of the ion are classified by each point group. Then, the 
non-interacting Hamiltonian Ho is given as 
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where is the annihilation operator of the ion with 
the symmetry index /i and energy and c^^ is anni- 
hilation operator of the conduction electrons with the 
radial wavenumber fc, angular momentum I = [l,m), 
and energy e^j. For the ion operators, the constraint 
0'l.o.f_t = 1 is required. 
The form of the interaction part of Hamiltonian _ffi,it 
is expressed in general as^ 
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where J~~ are coupling constants and "^7,,, are matri- 
ces of the 7-th irreducible representation (IR) of the 
direct product aj^cii, and subject to the constraint as 

y P = 1. 

For this kind of Hamiltonian H = Hq + i?i,it , 2-loop 
renormalization group (RG) equations have been derived 
formally in ref. 9. A set of RG equations is given as fol- 
lows: 



dx 
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where x = log(I?/£'o), D (Dq) being the half of the 
scaled (bare) conduction-electron bandwidth, and p is 
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Table I. Definition of 4?'^ = a^^'^^jj^^a,^ for (1, 0, 0) type cubic con- 
figurations, n^i = ajiafi, and '^m\r is a r component of p-th 
irreducible representation. 
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the density of states (DOS) at the Fermi level. For sim- 
plicity, wc take the same bandwidth and DOS for all the 
partial- wave components / of the conduction electrons. 
The RG equations, (4) and (5), are quite general and are 
not dependent on details of models. 

In Table 1, IR, , describing hopping between the 
different ion configurations for the six-level model, are 
listed. In six-level model, there are three IR's: r+(^ 
^t(-^ {ash^-si} = {0-3 (322-r2)/V3''*3 (x^-y^)}) and 
rj(-^ {a4x, a4y, a4z})- The energy of the IR's are es- 
timated as Ai < A4 < A3 for the case of retaining 
only the nearest neighbor hopping. For the conduction 
electrons, we retain the partial waves of / < 2 and de- 
fine S = J2k'^kOO,Pz = Z^feCfclo, V^Px = Z]fc(Cfel-l - 

Cfeii), V^py = J2k Kcfei-i+Cfcii), V^dxy = Ylk i(cfe2-2- 

C/c22), V^dyz = Z)fei(cA;21+Cfe2-l), V^dzx = Z)fe(cfe2-1- 

Cfc2i), e| = X^j. Cfc2o, and v^e^ = X^j. (0^22 +0^2-2)- The 
IR's for the conduction electrons -jmsr arc listed in Ta- 
ble 2. Using '^^j jmsr, we obtain the bare interaction 
Hamiltonian as 
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where we have used approximations a <C 1 and a local 
interaction parameterized by uq.^ The initial couplings 
J^/°^'s in the sense of RG are estimated from eq. (6). 

We solve the RG equations (4) and (5) numerically. 
The results of 1-loop and 2-loop RG are shown in Fig. 
1 and Fig. 2, respectively (we only show for typic;al cou- 
plings). The qualitative features are similar to the case 
of the four-level model. ^ In 1-loop calculation, we obtain 
the characteristic energy scale Tk/Dq ~ exp(— l/(/3'Uo))i 
at which some of the effective coupling constants diverge. 
In the 2-loop order, the coupling constants saturate at in- 
termediate values approaching the fixed point. What we 
have to be careful in Fig. 2 is that eqs. (4) and (5) are 
derived assuming D » max(Jj^, A,,). Then the reliable 
range of D/D„ is, at least, D/Da > A/Do - 0.001. It is 
noted that the energy of the first excited level of the ion 
decreases as renormalization proceeds, while that of the 
second one increases. 

In the case of four-level system, the correct fixed point 
corresponds to the results obtained by 1-loop calculation, 
which was verified by the NRG calculation. The detailed 
analysis shows that the s, Px and Py components of the 
conduction electrons form Kondo singlet with the ion 
with (singlet) and (doublet) symmetry imder D4/1 
point group. This can be seen in Fig. 3, in which we show 
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Table II. Definition of 7;i5;3r , r component of p-th irreducible rep- 
resentation Fl,. 
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Fig. 1. 1-loop RG flows. The parameters used are uq = 0.3-Do> 
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Fig. 2. 2-loop RG flows. The symbols and parameters are same 
as in Fig. 1. Inset: Effective levels of the ion states A^/Dq vs 
D/Do, M = 1,3 and 4. 



the phase shifts of the conduction electrons estimated by 
the NRG energy spectra. The orbitals of the s and p± 
strongly interact with the ion, while those of d± do not. 
These kinds of orbital Kondo effect seem to give an ori- 
gin of the heavy effective mass observed in SmOs4Sbi2. 
Indeed, we obtained the enhanced linear temperature co- 
efficient of the specific heat in this model by the NRG 
calculations.^ 

In the case of six-level system, the dominant coupling 
constants at low energy fixed point are those between 
(s, Px Py, Pz), and (r:;^, T^^, T^^, T^J, which is a 
natural extension of the results in the four-level sys- 
tem. At low temperatures, an orbital Kondo effect arising 
from these couplings are thought to become possible. Al- 
though we cannot carry out the NRG calculation for the 
six-level case because of the computational difficulty, we 
expect that the result of the 1-loop calculation captures 
the nature of the low energy fixed point as the four-level 
model. 

In the present paper, we have investigated (1, 0, 0) type 



Fig. 3. Phase shift of conduction electrons 5^ vs no. ^ indicates 
the orbital index of the conduction electron: p-|- = —{px + 
ipy)/V2, p- = (px - ipy)/V^, d+ = (e| + idxy)/V2 and 
d— = (e| — idxy)/ \/2. Lines are guide to eyes. 



six-level model. Recently, Kaneko et al. carried out neu- 
tron scattering experiments in PrOs4Sbi2, and deduced 
the nuclear density distributions.^" Their result suggests 
that the charge distribution of Pr ion extends mainly in 
the (1, 1, 1) direction at high temperatures and isotropi- 
cally at low temperatures. Concerning this point, it is de- 
sired to investigate the properties of an eight-level model, 
in which the stable points of the ion locate at eight (1,1,1) 
directions. We expect that the similar orbital Kondo ef- 
fects occur even in the eight-level model. 
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